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ABSTRACT 



We present predicted full-sky maps of submillimeter and microwave emission 
from the diffuse interstellar dust in the Galaxy. These maps are extrapolated 
from the 100/zm emission and 100/240/im flux ratio maps that Schlegel, 
Finkbeiner, & Davis (1998; SFD98) generated from IRAS and COBE/DIRBE 
data. Results are presented for a number of physically plausible emissivity 
models. 

The correlation of COBE/FIRAS data with the simple SFD98 (z/ 2 emissivity 
power law) extrapolation is much tighter than with other common dust 
templates such as Hi column density or lOOyum emission. Despite the apparent 
success of the SFD98 extrapolation, the assumed u 2 emissivity is inconsistent 
with the FIRAS data below 800 GHz. Indeed, no power law emissivity function 
fits the FIRAS data from 200 — 2100 GHz. In this paper we provide a formalism 
for a mult i- component model for the dust emission. A two-component model 
with a mixture of silicate and carbon-dominated grains (motivated by [Pollack 
et at. 1994] ) provides a fit to an accuracy of ~ 15% to all the FIRAS data over 
the entire high-latitude sky. Small systematic differences are found between the 
atomic and molecular phases of the ISM. 

COBE/DMR, has observed microwave emission that is correlated with 
thermal dust emission. However, this emission is higher than our model 
predicts by factors of 1.2, 2.4 and 20 at 90, 53 and 31 GHz, respectively. This 
provides evidence that another emission mechanism dominates dust emission at 
frequencies below ~ 60 GHz. 

Our predictions for the thermal (vibrational) emission from Galactic dust at 
v < 3000 GHz are available for general use. These full-sky predictions can be 
made at the DIRBE resolution of 40' or at the higher resolution of 6'.1 from the 
SFD98| DIRBE-corrected IRAS maps. 
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1. INTRODUCTION 

The pioneering Infrared Astronomy Satellite (IRAS) led to the discovery of the 
ubiquitous infrared cirrus, whose thermal emission is especially visible in the 100/zm band 
( [Low et al. 1981 ■ This cirrus, with a characteristic temperature of ~ 20 K, arches across 



the sky in long filamentary chains and is present at all Galactic latitudes. However, IRAS 
was optimized for the detection of point sources, and its ability to map the diffuse cirrus was 
less than optimal. Because of calibration drifts and hysteresis effects, the resulting IRAS 
Sky Survey Atlas (ISSA: [Wheelock et al. 1994Q images are contaminated by significant 
striping and poor control of large scale gradients. 

The Diffuse Infrared Background Experiment (DIRBE) on the COBE satellite is the 
perfect complement to IRAS. It has relatively low angular resolution (0.7°) but superbly 
controlled zero-points and gains. This has led to the generation of a map of the far-infrared 
sky with unprecedented accuracy and uniformity of coverage. Schlegel, Finkbeiner, & 
Davis (1998; hereafter |SFD98|) created a merged map of the IRAS and DIRBE data 



with an angular resolution of 6' and DIRBE-quality calibration. Their full-sky map 
shows the pervasive extent of the infrared cirrus and has proven successful for estimation 
of extragalactic reddening. But equally important will be the use of this type of data 
for estimation of Galactic foreground for the coming generation of CMBR experiments, 
including MAP and Planck and a host of ground- and balloon-based projects. 

In this paper, we consider the use of the SFD98 dust map as a predictor for 
microwave emission from Galactic dust. The SFD98 map is based solely upon 100 — 240/im 
(1250 — 3000 GHz) emission. Extrapolation to microwave frequencies is very sensitive to 
the details of the composition and emissivity properties of the dust. We show that the v 2 
emissivity assumed by SFD98 is inconsistent with the 100 — 2100 GHz emission probed 
by the COBE Far Infrared Absolute Spectrophotometer (FIRAS). We use these FIRAS 
data to constrain the properties of the dust, and show that no power law emissivity model 
can consistently explain the full spectral range of the dust emission. However, we find 
excellent agreement with a two-component model whose components we tentatively refer to 
as silicate and carbon-dominated grains. With this model for the dust emissivity function, 
extrapolation of Galactic dust emission from 100/im to lower frequencies is based upon the 
filtered DIRBE 100/240/im color temperature. 

In §^|, we discuss the COBE data sets and the details of comparisons using SFD98. 
Section § explores a variety of one-component dust models, demonstrating that a single 
power-law emissivity fails to explain the data, as does a broadened temperature distribution. 
Section |] explores a family of two-component dust models, in which energy balance and 
the temperature of the separate components are tightly coupled - one of which achieves 
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excellent agreement with the FIRAS data. Section [| discusses the robustness of this best 
model with respect to various ISM environments, and §[] compares our predictions to 
(DMR) microwave observations, demonstrating that the microwave emission may exceed 
the predictions of any thermal (vibrational) emission mechanisms. This is perhaps the 
signature of spinning dust grains emitting electric dipole radiation (|Draine fc Lazarian| 
1998b|) or the signature of free-free emission. Summary and conclusions are presented in §[7l 



2. DATA SETS 

The COBE (COsmic Background Explorer) satellite consisted of three 
instruments, DMR (Differential Microwave Radiometer), FIRAS (Far Infrared Absolute 
Spectrophotometer), and DIRBE (Diffuse Infrared Background Experiment). In this paper 
we shall compare predictions of dust emission based on DIRBE in the far-infrared with 
that observed by FIRAS at lower frequencies. In addition, we extend this correlation 
to still lower frequencies (31.5, 53, and 90 GHz) observed by DMR. Although the DMR 
fluctuations are dominated by intrinsic CMBR anisotropy, a residual correlation with 
DIRBE is detectable even at high latitudes. 



2.1. FIRAS Spectra 

The objective of the FIRAS instrument was to compare the cosmic microwave 
background radiation (CMBR) to an accurate blackbody, and to observe the dust and 
line emission from the Galaxy. It is a polarizing Michelson interferometer ( [Mather 1982; ), 
operated differentially with an internal reference blackbody and calibrated by an external 
blackbody with an emissivity known to better than 1 part in 10 4 . It covers the wavelength 
range from 0.1 to 10 mm (30 — 3000 GHz) in two spectral channels separated at ~ 0.5 mm 
(600 GHz). The spectral resolution is ~ 20 GHz. Although the design of the FIRAS 
experiment was optimized for its very successful measurement of the CMBR spectrum 
( [Fixsen et at. 1996 ), the instrument also measured the spectrum of the dust emission of our 



Galaxy (e.g., [Fixsen et at. 1996]). For the highest frequency channels, the Galactic signal 



dominates all others. 

A flared horn antenna aligned with the COBE spin axis gives the FIRAS a 7° field of 
view. The instrument was cooled to 1.5 K to reduce its thermal emission and enable the 
use of sensitive bolometric detectors. The FIRAS ceased to operate when the COBE supply 
of liquid helium was depleted on 21 September 1990, by which time it had surveyed the sky 
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1.6 times. 

We use the FIRAS Pass 4 Galactic dust spectra (hereafter FIRAS dust spectra) from 
which CMBR, zodiacal light, and a FIRB model have been subtracted (|Fixsen et al. 19971) . 
The data are presented as 213 spectral bins on the resolution 6 skycube map (6144 pixels 
on the full sky).[| Several Galactic emission lines, such as C + (157.7/xm), have been removed 
and replaced with interpolated values. For our analyses, we have removed the troublesome 
frequency bins listed in table Q, and recalibrated the entire FIRAS data set down by 1%. 
Our analyses make use of 123 frequency bins at 100 < v < 2100 GHz (140/xm < A < 3 mm). 
Note that data in the lowest two frequency bins are off the page in some of the figures, but 
are used in the fits. Details of the frequency bin choice and recalibration can be found in 
Appendix 

2.2. DMR Data 

DMR observed the sky at three frequencies, 31.5, 53, and 90 GHz, achieving the first 
detection of anisotropy in the CMBR ( [Smoot et al. 1992] ). In this paper we use the 4-year 
DMR skymaps dated 18 April 1995, which have the monopole and dipole removed. These 
maps do not influence any of our model fits, but are compared with our predictions in §|6]. 

Kogut et al. (1996) observed a correlation between Galactic dust and the 31.5 and 
53 GHz channels of DMR that is much greater than that expected from any models of 
thermal (vibrational) emission by dust. Alternative explanations such as spinning dust 
grains ( praine fc Lazarian 1998"fj| ) or spatially correlated free-free emission have been 
proposed, but are not well constrained by existing data (c./., |de Uliveira-Gosta et al. 19U8] ). 
We discuss this excess emission in 

2.3. DIRBE Data and SFD Dust Maps 

2.3.1. SFD emission map 

I I 

SFD98 presented a full-sky 100/im cirrus emission map constructed from both the 
DIRBE and IRAS/ISSA data sets. The map is well calibrated, zodiacal light subtracted, 
Fourier destriped, and point source subtracted, with a final resolution of 6.1'. Complete 



1 These data are available on the World Wide Web at http://www.gsfc.nasa.gov/astro/cobe/. 
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descriptions of these maps may be found in |SFD98 .^ For comparisons with FIRAS and 
DMR, the full resolution of the IRAS/DIKBE map is not required. Instead, we use the 
0.7° DIRBE map with zodiacal light removed as described in SFD98, with point sources 
included. The DIRBE map offers a fair comparison with the high-frequency FIRAS data, 
in which these sources contribute to the measured flux. The comparison is less appropriate 
in the low-frequency FIRAS data, where typical FIR-luminous sources are expected to 
contribute little to the measured millimeter flux. However, the contribution from stars and 
galaxies to the 100/im flux is only ~ 2% of the diffuse Galactic emission at high latitudes, 
and relatively less at low latitudes. 



2.3.2. Ratio Map 

We also make use of a DIRBE 100/!m/240/im color ratio similar to that described in 
SFD98. Because of the poor signal-to-noise ratio in the 240/im map, SFD98 employed a 
filtering algorithm to give the minimum variance estimate of the dust temperature in each 
1.3° Gaussian beam. In each pixel this filter yields the weighted average of the measured 
flux and a more robust estimator - in the case of SFD98, the estimator is the |6| > 75° 
average flux. The weights are chosen so that the ratio of the filtered maps is the minimum 
variance estimate of the true flux ratio. The process gives the measured ratio in high S/N 
pixels, but recovers the high latitude average ratio of 0.66 in the limit of low S/N. 

The SFD98 algorithm has the unfortunate effect of suppressing temperature variations 
at high latitude even when those variations are measurable at a resolution of a few degrees. 
In the current analysis, the S/N of FIRAS in a 7° beam is sufficiently high that this 
non-local filtering algorithm causes undesirable behavior in model fits. When the DIRBE 
100/im and 240/im maps are smoothed to 7°, structure in the ratio appears which is not 
aligned with the imperfectly-subtracted zodiacal plane or other potential artifacts in the 
maps. Rather, the DIRBE 240/im map exhibits structure, even at very low levels, that is 
correlated with the FIRAS maps at 240/im. Therefore, it is presumed that this structure is 
of extra-solar origin, and should not be discarded as it was in the SFD98 analysis. 

In the current paper, we have constructed a new ratio map, TZ, that retains more 
temperature information. We use the same weight function W described in SFD98 
(equations 8 and 9). But rather than forcing the map to a high-latitude average at low S/N, 
we force it to the local 7° average. High S/N regions are little changed from the previous 
R map, but large-scale temperature structures are now apparent at high Galactic latitudes 



2 These data are publicly available via the World Wide Web at http:/ /astro. berkeley.edu/dust. 
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that were suppressed before. The temperature correction derived from this ratio map have 
a 1.3° resolution in high S/N regions, and is applied to the full- resolution 100/i map, not 
the smoothed 100/im map. This procedure correctly handles the situation where a compact, 
high S/N source is located near a diffuse background with a different color temperature. It 
should be noted here that the same 13 bright sources listed in SFD98 Table 1 were removed 
from the DIRBE maps before smoothing, to avoid halo artifacts in the 1Z map. This change 
in 1Z produces only very modest change in the SFD98 reddening predictions. The largest 
change to predicted reddenings at high latitude is of order E(B — V) = 0.01 mag. 



2.3.3. Cosmic IR Background Removal 

The ratio map and derived temperatures are moderately dependent upon the 
uncertainties in cosmic infrared background (CIB). The CIB represents the extragalactic 
signal that is unresolved and isotropic in either the DIRBE or FIRAS instruments. This 
signal is presumably from high-redshift [z ^ 1) dust-enshrouded galaxies, which are only 
beginning to be resolved with ground-based submillimeter observations (c.f., [Blain et aU 
1999| ). Detections of the CIB at 140 and 240/im were reported last year by |SFD98| and 



Hauser et al. 199"E| , and upper limits were reported at 100/im. A more definitive analysis 



is in preparation by Finkbeiner, Schlegel, & Davis (1999). We remove the CIB from the 



DIRBE maps in the same way as (5FD98| - as part of the zodiacal light model. By using the 
zero point of the Leiden-Dwingeloo Hi map, a model including CIB and zodiacal light may 
be fit and removed. This is an easier problem than the separation of CIB from zodiacal 
light, which is unnecessary for this paper. One source of error in this could result if there is 
significant dust emission correlated with Ha, and the Ha/dust correlation has a different 
zero point than the Hi /dust correlation. The sense of this would be to add a constant 
to both Jioo and J240, causing the derived temperature distribution on the sky to broaden 
or narrow. In other words, a poor TZ would produce FIRAS fit residuals that depend on 
temperature, which would show up in figure |^. Lack of a temperature-dependent residual 
indicates that CIB removal errors have had a negligible effect on our model. The rest of our 
fit procedure - using correlation slopes at each frequency - ignores an isotropic component 
by construction, so we conclude that we are unaffected by uncertainty in the CIB. 



2.4. Comparing COBE Data Sets 
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2.4-1- Beam Shapes 

Comparisons between DIRBE and FIRAS data are made at the FIRAS resolution. 
The FIRAS beam has a frequency- dependent shape that is not well-measured. The beam 
is approximately a 7°-diameter tophat in the highest-frequency channels, with power-law 



wings (measured at 750 GHz from off-axis measurements of the moon; see §7.9.4 of [FIRAS 
[Exp. Supp. 19971) . The beam shape is closer to Gaussian at lower frequencies, with 
exponential wings from 5° to 15° from the beam center (measured at 90 GHz in the lab; 
see Fixsen et al. 1994] ) . Because the FIRAS scan strategy averages over 32 to 46 sec 
interferograms, the beam is smeared by typically 2.3° in approximately lines of constant 
ecliptic longitude. The pixelization of the FIRAS data on ~ 3° pixels introduces another 
effective smoothing. We approximately match the FIRAS beam by first convolving the 
DIRBE data with a 7.0° circular tophat, then convolving with a 3.0° circular tophat, then 
smoothing by 2.3° in ecliptic longitude. We ignore the non-Gaussian beam shapes of the 
DIRBE instruments since they are sufficiently smaller. 

We attempted to match the frequency-dependence of the FIRAS beam. The signal 
in the Galactic plane is sufficiently strong (~ 100 times larger than the median value 
at 500 GHz) that the exact sidelobe profile may be important. At low frequencies, the 
sidelobes exceed 10~ 3 within 8° of the beam center. However, since the profile has only 
been measured at two frequencies, it is impossible to model the beam to high accuracy. 
Therefore, the beam shape uncertainties introduce errors of up to 10% within 7° of the 
plane. Because we exclude the sky within 7° of the Galactic plane from our analyses for 



other reasons (see § [2.4.2|) , we simply ignore the complication of frequency-dependence of 
the FIRAS beam. 



2-4-2. Spatial Mask 

Our analysis is limited to those parts of the sky where the far-infrared emission is 
expected to be dominated by the diffuse interstellar medium. We create a spatial mask that 
excludes the Galactic plane below \b\ = 7°, the Magellanic Clouds, and Hn regions in Orion 
and Ophiuchus. In such regions, the SFD temperature map is unreliable due to confusion 
limits. These are also the regions where the FIRAS data suffer from poorly understood 
sidelobe contamination. We also mask 1.3% of the sky where the FIRAS coverage is missing 
or incomplete, and another 15% where the FIRAS pixel weight is less than 0.4 (the median 
value is 0.8). The final mask excludes 29% of the sky from our analyses, and is shown as 
the thin black outlines in Figs. ||, |9|, and This mask is used throughout this paper 
except for the comparison in figure p] in which the Galactic plane is included, and for the 
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comparisons with DMR shown in Tables |3] and |]. For the DMR comparisons we apply the 
Goddard "custom cut" mask from the 4-year DMR data analysis which excludes 37% of 
the sky QBennett et at. 19961 ). 

2.4-3. Simple Difference Spectra 

For an overview of the three COBE datasets, it is useful to plot the DMR, FIRAS, and 
low-frequency DIRBE data on a single plot. The DMR is a differential instrument, so the 
mean measurement over the sky is zero in each channel. In order to compare it to FIRAS 
and DIRBE, we plot the difference between "bright" (/ 900 gh z > 3.0 MJy/sr) and "faint" 
regions of the sky. This method has the additional advantage of discarding any isotropic 
background of cosmic or instrumental origin. We have further divided the sky into cold, 
warm, and hot zones according to the DIRBE iioo/^o color ratio. The cold component 
(/ioo/^24o < 0.62) comprises 14% of the sky, the hot component (/ioo/-^24o > 0.69) comprises 
26%, and the warm component comprises 44%. The remaining 16% is masked, rejecting 
only bad or noisy FIRAS pixels, but including signal in the Galactic plane for better S/N. 
This plot assumes a monotonic relationship between color ratio and physical temperature, 
but requires no other knowledge of the dust spectra. The difference spectra for these regions 
are shown in Figure |l|. All three curves are renormalized such that the DIRBE lOOyum 
flux is 1 MJy/sr, a value typical for high-latitude dust. The spectra qualitatively have 
the correct behavior, with the "cold" regions showing stronger emission at low frequencies 
relative to 100/im. 

The FIRAS emission at low frequencies (200 <^ v <^ 600 GHz) scales as ~ u 3 ' 2 . Because 
the Planck function B V (T) asymptotes slowly to v 2 on the Rayleigh- Jeans tail, the best fit 
emission spectrum is ~ u l 7 B u (19 K) over this frequency range, not ~ i/ 1,2 -B !/ (19 K). The 
temperature 19 K corresponds to the median dust temperature for this fit to the emission 
spectrum. These considerations indicate that the naive v 2 emissivity law assumed in SFD98 
is incorrect, a matter that will be explored extensively in this paper. 

Furthermore, Figure [I] demonstrates that the Galactic emission detected by DMR is 
inconsistent with any power-law extrapolation of the FIRAS data. The DMR 31.5 and 
53 GHz channels lie well above the power law extrapolation of the FIRAS curves. We 
address this problem in 



3. ONE-COMPONENT DUST MODELS 
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3.1. Predicted Microwave Emission from SFD98 

A simple but naive prediction for submillimeter /microwave emission can be made 
from our previous work. [SFD98| extensively studied the emission from dust in the regime 
100/im < A < 240/im. Assuming a v 2 emissivity model, the temperature of the dust was 
mapped with a resolution of 1.3° from the DIRBE 100/im/240/im emission ratio, TZ. The 
100/im emission of the dust was mapped with a resolution of 6.1' by utilizing small-scale 
information from the IRAS mission. Emission at lower (millimeter/microwave) frequencies 
can be predicted by extrapolating the 100/im flux using this temperature fit. For each 
line-of-sight in the maps, the emission at frequency v can be expressed as 

L = K w l(a,T)I 1Q0 ^^ (1) 

where B V (T) is the Planck function at temperature T, iioo is the DIRBE-calibrated 100/im 
map, Kioo(a,T) is the color correction factor for the DIRBE 100/im filter when observing a 
u a B u (T) spectrum, and u = 3000 GHz is the reference frequency corresponding to 100/im. 
Our values for the color correction factor can be recovered for all values of a used in this 
paper from the formula 

where the a and b coefficients may be found in Table [|. 

The choice of an a = 2 emissivity model was not well-motivated in [SFD98|. The dust 



column map is only very weakly dependent upon the emissivity law because the entire 
map is renormalized using direct observations of reddening. Using an a = 1 emissivity 
model changes the relative column density of dust between warm and cold regions by only 
~ 1%. However, the extrapolated emission at lower frequencies is highly dependent upon 
the emissivity of the dust, and the a = 2 assumption must be tested. 



3.2. Comparison with FIRAS Data at 500 GHz 

The extrapolated, millimeter emission from dust as predicted from |SFD98| can be 
compared directly to the FIRAS measurements. As a test of this model, we examine the 
spatial correlation of FIRAS 500 GHz emission with the DIRBE 100/im map. Such a 
comparison is most meaningful at a FIRAS frequency that is nearly on the Rayleigh- Jeans 
tail of the dust spectrum, but is still easily measured against the 2.73K CMBR. This 
strikes a balance between the poor S/N ratio at lower frequencies and the stronger 
temperature-dependence of the spectrum at higher frequencies. 
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To increase S/N, we synthesize a broadband FIRAS 500 GHz channel. We sum the 
FIRAS measurements in the 10 unmasked channels, i, between 400 and 600 GHz, weighting 
by z/~ 3 - 5 to make the summand roughly constant at these frequencies: 



J2i /•>, '■ 



FIRAS 500 = ^ ; J3.5 ( 3 ) 



The (5FD98| extrapolation (equation |T|) is similarly summed over the same FIRAS frequency 



bins to generate a predicted broadband flux: 

SFD 500 = o » ' ol j 4 

The correlation between FIRAS 500 and SFD 500 is very good (Figure |2|c) , with an RMS 
dispersion of 0.2 MJy/sr about the best-fit line. For comparison, we plot the correlation 
with (a) Hi column density, and (b) DIRBE 100/im flux. Both of these show a scatter that 
is 3.7 times worse than with SFD 500, demonstrating that sub millimeter emission from 
dust is neither simply related to the Hi column density, nor is the dust at one temperature 
everywhere on the sky. 

The SFD98 extrapolations work impressively well in predicting 500 GHz emission from 
dust, despite their assumptions of v 2 emissivity and one temperature along each line-of-sight. 
However, the slope of the regression between FIRAS 500 and our extrapolations differs 
significantly from unity (formally by nearly 40a). At lower frequencies, the slope departs 
even more strongly from unity. This is an indication that a v 2 emissivity is incorrect for the 
dust, as was seen from the mean spectrum of large regions of the sky (Figure [I]). This will 
be addressed in detail in 5B.3L 



3.3. The Spectrum of Dust-Correlated Emission 

The many frequency channels of the FIRAS experiment allow detailed comparisons 
with predictions for the spectrum of dust emission. For each channel of the FIRAS data, 
we compute a correlation slope with the SFD prediction. The correlation slope is computed 
as the best-fit slope of the FIRAS column (F p ) versus the predicted column (I p ). By 
subtracting a weighted mean from each map, the correlation slope, m, is insensitive to 
zero-point uncertainties in either map: 

(W P (F P -F P )(I P -I P )) 

(W P (I P - I p ) 2 ) 1 ; 
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where W p is the FIRAS pixel weight (W p ~ 1/^) for pixel p. Such a slope is computed for 
each FIRAS channel centered at v^. These slopes are equivalent to that computed in Figure 
|2]c for a broad-band FIRAS channel. The correlation between the FIRAS column and the 
SFD prediction is strong and apparently free from systematic errors in all but the lowest 
frequency channels. For all FIRAS channels v < 2100 GHz, we compute m using 71% of 



the sky described in §2.4.2. If the v emissivity model used by SFD98 were valid, the slope 



would be consistent with unity. 

In each panel of Figure 0, the correlation slopes, or "slope spectrum," are plotted as 
a function of frequency, with the vertical lines extending to ±3cr. Overplotted in each 
panel are various models, evaluated for a typical high-latitude 1Z value of 0.68. To facilitate 
comparison with the data, these models are each divided by the same v 2 prediction as the 
data. Such a comparison is instructive, but not used for formal fits, because the ratio of a 
given model prediction to the v 2 prediction depends weakly upon 1Z, and therefore varies 
across the sky. Note that all model ratios are constrained via our temperature correction 
to be unity at 100 and 240/mi (3000 and 1250 GHz). The data points are also constrained 
to go through unity at 1250 GHz to the extent that the FIRAS and DIRBE data are 
consistent. 

At v > 500 GHz, this slope spectrum is consistent with the v 2 model to within 10%. At 
lower frequencies, the slope spectrum increases, demonstrating that there is more emission 
at low frequencies than a v 2 emission model would predict. 



3.4. Other Power-Law Emissivities 

There is no power-law emissivity model that fits the FIRAS data. The SFD prediction 
can be made with other emissivity profiles by modifying the exponent, a, in equation [l]. 
An a — 1.5 emissivity profile results in a better fit at low frequencies, but ruins the fit at 
high frequencies (see Figure ^b). An a = 2.2 emissivity gives a good fit at high frequency, 
but is catastrophically wrong at low frequency. The minimum y 2 is achieved for an a ~ 2.0 
extrapolation. The value of x 2 doubles for a = 1.5 (see Table [I]). 



3.5. Broadened Temperature Distribution 

Our model ignores the possibility of dust temperature variation along a line-of-sight 
through the Galaxy. Such a situation may arise from the superposition of different 
environments with different temperatures. Or it may arise from an intrinsic distribution in 
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dust grain sizes and temperatures within a given environment. 

The far-IR/submillimeter dust emission is expected to be dominated by large grains 
(O.Olyum < a < 0.25/im), which are in equilibrium with the ambient radiation field. The 
grains are thought to follow a power law distribution of grain sizes dn/da oc a~ 3 ' 5 ( [Mathii 



et al. 19771) from about 0.0005/mi (the size at which absorption of a single photon can 
sublimate mass away from the grain) to 0.25/im, where the number density appears to fall 
off based on Ry measurements (|Kim, Martin fc Hendry 1994]) . Since the long wavelength 



emissivity of a grain scales as its size a times the surface area, or volume a 3 , the larger grains 
dominate the submillimeter emission ( praine fc Anderson 1985|) . This will be true unless the 
slope of the size distribution is steepened to a slope of nearly —4.0. The very small grains 
(VSGs; a O.Olyum) are transiently heated and emit at high effective temperatures for a 
small fraction of the time, but do not contribute significantly to submillimeter emission. 

Even for the grains which are large enough to be in equilibrium with the ambient 
radiation field, there is a slight size-dependent temperature variation. The approximation 
that the grains are small compared to the wavelength of absorbed radiation is not exactly 
satisfied at the large end of the grain size distribution, so the larger grains are a bit colder 
because they absorb less efficiently relative to their emission. For reasonable assumptions 
about the ISRF, the temperature varies approximately as a -0 06 ( praine fc Lee 1984| ), both 



for silicate and graphite grains. Over the size range of interest, O.Ol/i < a < 0.25fi, the 
temperature range of the grains at a given locale is modest. The dominance of the largest 
grains' emission results in a narrow distribution of relevant grain temperatures, and allows 
us to use the emission- weighted mean temperature for each component. This approximation 
is good at 100/im and exact on the Ray leigh- Jeans tail where I u oc T. This greatly simplifies 
our analysis. 

There are two other reasons why the temperature might vary along a line of sight. The 
cutoff of the grain size distribution at large a might vary, causing the dominant size (and 
temperature) component to vary along the line of sight. Another possibility is that the 
interstellar radiation field may vary. Along lines of sight passing through cold molecular 
clouds, both of these effects should contribute. For extinction predictions (as in SFD98) it 
is important to understand which of these effects is causing the temperature variation; for 
extrapolating the emission to microwaves, the cause is unimportant. 

In order to model such variations, we experimented with Gaussian-broadened 
distributions of temperatures with width AT FWHM along a single line of sight. The 
ratio of a broadened-T model (ATfwhm = 3 K) to the idealized single-temperature fit is 
shown as the dashed line in Figure |3|a. An even broader distribution (ATfwhm = 6 K) is 
plotted as a dash-dot line. The broadened-T model changes the predictions by at most 
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4% for AT FWHM = 3 K and 15% for AT FWHM = 6 K. These models retain large and 
coherent deviations from the FIRAS data. The models are slightly more consistent with 
the low-frequency FIRAS channels, but less consistent with the high-frequency channels. 
The value of x 2 f° r these models is higher than that of the single-temperature v 2 model. 
We conclude that a Gaussian-broadened temperature distribution does not fit the data 
any better than a single-temperature model and would only introduce poorly-constrained 
parameters into our models. The lack of an acceptable one-component model - even 
when temperature variation along each line of sight is included - indicates the need for 
multi-component models, which are discussed in the next section. 



4. MULTI-COMPONENT DUST MODELS 

In this section we explain the theoretical motivation for multi-component dust models, 
present our general model, and then provide the results for eight specific models. 



4.1. Theoretical Motivation 

The diffuse ISM is known to contain many different types of molecules and dust grains 
with a broad range of physical properties. In spite of the expected melange of dust grains, 
it was originally expected that in the far-IR/sub millimeter bands, all dust would have 
similar optical properties. For example, Draine & Lee (1984) predicted v 2 emissivity for 
both silicate and graphite grains. 

The emission mechanism corresponding to fundamental vibrations (single 
photon/phonon interactions) in crystalline dielectric materials is optically inactive 
due to wavevector conservation, and multi-photon interactions are rare at low temperatures. 
Therefore, the emissivity of crystalline materials would be expected to be dominated by 
absorption in the damping wing of an infrared active fundamental vibration, the strength of 
which goes as v 2 at low frequencies. In a metallic or semimetallic material, interaction with 
electrons was expected to dominate FIR absorption, also resulting in v 2 emissivity ( |Wooten 
19720 . 

In amorphous materials, the lack of long range order causes a breakdown of the 
selection rules that forbid single photon/phonon interactions, and all modes become active. 
The emissivity power law then only depends upon the density of states, which was also 
thought to go as v 2 ( [Kittel 1976|) . Thus, amorphous materials were expected to have the 
same dependence on frequency as other components, but for an entirely different reason. 
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The most notable exception to the v 2 theory was the case of planar structures such as 
graphite, which would yield a v 1 power law by the same reasoning. For an excellent 
summary of the theoretical details, see [Helens fc Allamandola (1987)| . 



More recent laboratory measurements suggest that universality of v 2 emissivity is 
an oversimplification, with different species of grains having differing emissivity laws. 
The composition and abundance of grains of different species can be constrained by 
astronomical observations and by observations of solar system bodies. A multi-component 
model for interstellar dust has been constructed by Pollack et al. (1994), based on 
laboratory measurements, observations of molecular cloud cores, and fits to dust shells with 
temperatures T rs 100 K around young stars. Their model predicts that at frequencies 
v ?Z 500 GHz, dust emission will obey a v 2 ^ emissivity law due to the dominance of carbon 
species. At lower frequencies, the emission is dominated by astronomical silicates such as 
olivine ([Mg,Fe]2Si04) and orthopyroxene ([Mg,Fe]SiOs). This low-frequency prominence of 
silicates flattens the emissivity profile to z/ L5 at frequencies v ^ 500 GHz (A <^ 600/im) 

Despite the complexities in the dust composition, most authors have chosen for 
simplicity to model the observed emission with a single power-law emissivity. If one's 
observations are limited to less than a decade in frequency, this parameterization may be 
adequate to fit the data, especially if one component dominates the emission. However, 
combining data from all three COBE instruments results in a tremendous range in observed 
frequencies. The discussion in demonstrates that a single power-law emissivity is a 
poor fit to this combined data. Our physical interpretation is that different grain species 
dominate the emission at different frequencies. 

Laboratory measurements of submillimeter-wave absorption properties of both 
crystalline and amorphous silicates ( |Agladze et al. 1996| ) suggest that a ranges from 



approximately 1.2 — 2.7, with some components having a much higher opacity than others 
(by a factor of ~ 40 at 300 GHz and 20K). These studies motivate a broader search of 
parameter space. 



4.2. General Multi-Component Model 

We outline a general formalism for describing a mixed population of dust grains. These 
simple considerations apply only in the limit of large grains (a > 0.01/im) which are not 
transiently heated but instead reach equilibrium with the local radiation field. We neglect 
emission from very small, transiently heated grains because it is unimportant over the 
FIRAS frequency range. 
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4-2.1. Statement of Model 

Fortunately, for wavelengths A > 100/x the large grains totally dominate the thermal 
(vibrational) dust emission. Since it is common to assume that each component of the dust 
will have a power-law emissivity over the FIRAS frequency range ( Pollack et al. 1994 ), one 



can sum these to construct a multi-component model analogous to equation |l|: 

T J2k fkQk{v)B u (T pk ) 



J2k fkQk{vo)B vo (T pk )K wo (ak, T pk ) 



where f k is a normalization factor for the k-th component, T pk is the temperature in pixel p 
of component k, K wo is the DIRBE color-correction factor ( PIRBE Exp. Supp. 1995|) and 
Ip,wo is the SFD98 100/mi flux at pixel p in the DIRBE filter. The emission efficiency Q{v) 
is the ratio of the emission cross section to the classical cross section of the grain. Because 
the grains of interest are very small compared to the wavelength of emission, Q(v) oc a, 
where a is the radius of a spherical grain. One interpretation of this is that the grain is so 
small that all parts of the grain are close enough to the grain surface to take part in the 
emission. 

The emission opacity (effective area per mass) for a spherical grain of radius a, K em [y\ 
is related to Q{y) by 



Because Q/a is usually taken to be independent of a for a « A (c.f., [Hildebrand 1983| ), 



K em (u) does not depend on grain size. The frequency dependence is taken to be a power law 

K em (u)=K em (u )(u/u )«" (8) 

where a k is the emissivity index and K, em (i> Q ) is the opacity of species k at a reference 
frequency uq = 3000 GHz. It will be convenient to interpret f k as the fraction of power 
absorbed and re-emitted by component k, so we force the power fractions to sum to unity: 

= (9) 

k 



4-2.2. Temperature Coupling 

The degrees of freedom in this multi-component model can be substantially reduced 
by demanding that the components are in equilibrium with the interstellar radiation field 
(ISRF). If we assume that the ISRF has a constant spectrum everywhere on the sky and 
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varies only in intensity, then we may define K* k to be the effective absorption opacity 
(cross-section per mass) to the ISRF in the limit of low optical depth: 

JO IlSRF{V)dv 

where IisrfW) represents the angle-averaged intensity in the ISRF as a function of 
frequency, and has the same dimensions as B v . To avoid confusion, k^ s (u) designates the 
optical opacity of component k, which is physically related to the submillimeter opacity 
K| m (V) but need not be an extension of the power law expression for /t| m (z/) in equation |8|. 
The total power absorbed per mass for species k is given by 



oo 



Ul n = 4 IisrfWv. (11) 







The power is primarily radiated in the far-infrared, and thus is only sensitive to the far-IR 
emissivity law. The power per mass emitted is 

U° k ut = / KT{v)B v {T k )dv. (12) 
Jo 

Demanding that each grain species is in equilibrium with the ISRF (e.g., Uj* 1 = U£ ut ), the 
energies of all species are related via: 

1 roo 1 poo 

- Kr(v)B v (T % )dv = - / K™(v)B v (T 3 )dv. (13) 

K* JO Kj Jo J 

Using our parameterization of the emissivities (equation^), we can solve for the temperature 
of one component clS db function of the other: 



qiZ(a>i) \ k B J 3 

where 

q = K em (v )/K* (15) 

is essentially the ratio of far-IR emission cross section to the UV/optical absorption cross 
section, and the integrals are absorbed into the analytic function 

Z(a)= -dx = C(4 + a)r(4 + a). (16) 



Henceforth, we shall use only the ratio of opacities, q k , assuming that the dust temperature 
is sensitive only to this ratio of emission to absorption cross sections. This is not strictly 
true, because k* is weakly dependent upon grain size. However, as we showed in § [375] , the 
assumption of a single temperature for each component - and therefore a single q% - in each 
locale is justified. 
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4-2.3. Interpretation of f , q 

Each dust component is therefore described by three global parameters (fa, q^, and 
one parameter that varies with position on the sky, T k (x) . Because equation [TJ| couples the 
temperature of each component, there is only one independent temperature (e.g., T2) per 
line-of-sight. The interpretation of the q k as IR/optical opacity ratios is obvious, but the 
meaning of the normalization factors is less clear. To understand what the /& are, let us 
consider the ratio of the power, Pi/Pj, absorbed and re-emitted by components i and j: 

P _ fi f Q °°^/^B u (T t )du 

and [15|, we see that the integrals are equal: 

Pi/Pj = fi/fj- (18) 

Because the fa sum to unity (equation |), we identify fa with the fraction of power absorbed 
from the ISRF and emitted in the FIR by component k. Note that Pi/Pj is independent 
of frequency. Note also that fa/^t * s proportional to the mass fraction. Therefore, if the 
optical opacities of all species were equal (which is unlikely), then fi/fj would measure the 
mass ratios between species i and j. 

Whether or not the actual components of the dust physically correspond to these 
components, equation || can be thought of as a phenomenological "expansion set" for 
describing the composite dust spectrum. 



Combining equations BL 13 



4.3. Fit Results for Specific Models 

We now describe our fits to eight different models of the form described above. All 
results from these fits are described in Table |I[ In the first four we consider only a single 
component (f\ = 1) whose temperature varies on the sky. We strongly emphasize that none 
of our models force a constant temperature everywhere on the sky, which is an extremely 
poor description of the data (as can be seen in Figure ^p). For a given a dust model, the 
temperature is uniquely constrained by the DIRBE 100/zm/240/im flux ratio along each 
line of sight. We perform single-component fits for a = 1.5, 1.7,2.0, and 2.2, obtaining the 
best reduced at a = 2.0. This is in agreement with previous fits to the FIRAS data 
( |Boulanger et al. 1996|) and also with earlier theoretical prejudice ( Praine fc Lee 1984j) . 



These results are encouraging, but statistically af = 3801 for 123 degrees of freedom 
(for the v 2 model) is completely unacceptable. The first two-component model we consider 
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is one designed to replicate the spectrum in Pollack et al. (1994). The prescription in their 
paper for their best-fit broken power law with (Ti) = (T 2 ) corresponds in our model to 
a\ = 1.5, a 2 = 2.6, f\ = 0.25, and qx/q-z = 0.61. This results in a considerably better fit of 
xl = 15.3 without fitting any new free parameters. The choices for ax, ct 2 , fx and qi/q2 
are based upon other empirical evidence completely independent of the DIRBE and FIRAS 
data sets. The ratio of the Pollack et al. model to the strawman v 2 model is shown in 
Figure |3]c as a light solid line. Between 800 and 1800 GHz, where the FIRAS signal is very 
good, the model matches the data to approximately 1% everywhere. At lower frequencies, 
the largest deviation is 25%. 

Allowing fx and qx/q2 to float with fixed «i,a 2 provides even better fits. We attempt 
to reproduce the results found in Reach et al. (1995), where a component of very cold dust 
was proposed to explain the low-frequency excess. Reach et al. used a = 2 for both the 
warm and cold components as a mathematical convenience. Letting our model float with 
«i = «2 = 2 we obtain fx = 0.00261 and 51/(72 = 2480. This model achieves a better fit 
than any other model tested in the literature do date, yielding a \ 2 = 1241, or 10.3/DOF. 
A physical interpretation of this combination of parameters in the context of our models 
would imply that there is a component constituting 0.26% of the dust emission power, but 
with an opacity ratio 2480 times higher than the dominant component. This huge opacity 
ratio explains the low temperature, (Ti) = 5 K, as compared to the dominant component 
at (T 2 ) = 18 K. Since this compelling fit appeared in the Reach et al. (1995) paper, some 
authors have sought to explain the model in the context of this simple interpretation. 



Fractal grains QFogel fc Leung 1998|) and other possibilities have been raised to explain the 
opacity ratio, but no convincing mechanism has yet been proposed. Opacities may indeed 
differ considerably, but factors of many thousand are probably unreasonable. However, 
the idea of multiple well-mixed components at different temperatures deserves further 
exploration. 

By taking the values from Pollack et al. (1994) of ax = 1.5, a 2 = 2.6, but realizing 
that such different components are very likely to be at different temperatures, we allow 
fx and qx/q 2 to float, obtaining X 2 = 244 or 2.0/DOF for f\ = 0.0309 and qx/q 2 = 11.2. 
It may seem surprising at first that one component is 11 times "better" at thermally 
radiating than another, but to justify this we appeal to the empirical measurements of 
Agladze et al. (1996). They find that the amorphous silicate MgO-2Si0 2 at 20 K, for 
example, radiates ~ 40 times more readily at 300 GHz than the crystalline silicate enstatite 
(MgSiOa) (|Agladze et al. 1996|, Table 1). The effective optical opacities k* may also vary 



significantly, so a wide range in emissivity ratios g 2 /gi is empirically well established. 
Furthermore, our model only requires a tiny fraction of the dust to be of this kind. This is 
a very reasonable theoretical step to take to obtain a formal increase in likelihood of ~ 350 



-19- 



orders of magnitude over a simple v 2 model. 

A further reduction of x 2 — 219 or 1.85/DOF is achieved by allowing the power law 
indices a.\ and a 2 to vary. The best-fit values, ati = 1.67, a 2 = 2.70, fi = 0.0363, and 
q 1 /q 2 = 13.0. 

For these model parameters, the temperatures of the two components are related by 

Ti = 0.352 T 2 118 . (19) 

The mean temperatures are (Ti) = 9.4 and (T 2 ) = 16.2 for the 71% of sky that we fit. This 
is the model we adopt for the comparisons discussed in the next section. 



5. DISCUSSION 

5.1. Interpretation of Best-Fit Model 

The thermal emission from Galactic dust can be very successfully predicted at 
millimeter /microwave frequencies using a two-component composition model with 
temperature varying on the sky. We tentatively refer to the two components as an 
amorphous silicate-like component [y x l emissivity, (T) pa 9.5 K) and a carbonaceous 
component (z/ 2,7 emissivity, (T) pa 16 K). This solution agrees with the FIRAS data much 
better than a 2-component model using two z/ 2 emissivity components with one of the 
components very cold ((T) ~ 5 K) Note also that we have obtained only 4 global model 
parameters from the entire FIRAS data set - all other column density and temperature 
information is derived from the DIRBE 100/im and 240/zm maps. In the Reach et al. 
analysis, temperatures of the two components were allowed to float independently, and were 
fit directly to the FIRAS data. 

Although our analysis does nothing to rule it out, we find no evidence for a recently 
proposed warm component (z/ 1 emissivity, (T) pa 29 K) associated with the WIM. This 
component results from a different approach to modeling the dust emission spectrum (see 
Lagache et al. 1999 for details). 



5.2. Spatial Coherence of Dust Properties 

Previous two-component models of dust emissivity fit to the FIRAS data found that 
the two components must be spatially correlated to a high degree ( [Reach et al. 1995]) . We 
demonstrate this fact by computing the correlation slope of dust with FIRAS 500 GHz as 
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shown in Figure |2|. Removing this correlated emission reduces the variance in the 500 GHz 
map by 95%. Based on the data at high Galactic latitude, about 4% of the variance is 
attributable to measurement noise. This leaves 1% of the variance (or 10% of the signal) 
as the upper limit for uncorrelated 500 GHz emission. If there exists a separate cold dust 
component which does not emit at 100/im, then it must be very highly spatially correlated 
with warm dust. 



5.3. Evidence for Variations in Dust Properties 

The agreement between the best-fit two-component model and the FIRAS data is 
impressive (see Figures | and ^). The reduced xt is close to unity, implying that the model 
uncertainties are small compared to the measurement errors. The far-IR/millimeter sky at 
\b\ ^ 15° appears to be well-fit by a fixed model for the interstellar dust. The temperature 
varies, but the composition and size distribution of the dust grains are constrained to be 
very similar everywhere in the diffuse ISM. 

Although our best-fit model appears to successfully describe the average dust emission 
spectrum, there still might be systematic variations across the sky. Splitting the sky into 
different zones based upon various observables, one can search for regions that deviate from 
our model. Dividing the sky according to temperature or dust column density, we do not 
find significant differences (see Figure |5|). 

However, splitting the sky according to dust/gas ratio, one does find coherent 
differences. We construct a dust/gas ratio by dividing the SFD98 dust map by the 
Leiden-Dwingeloo Hi map ( [Hartmann fc Burton 199"B] ), both smoothed to a 45' FWHM 
Gaussian beam. The SFD98 dust map is proportional to the 100/im emission expected if 
the entire sky were at a uniform temperature. Regions where the dust/gas ratio exceeds the 
high-latitude average by more than a factor of 2 are designated "molecular" (14% of the 
sky). Remaining pixels are designated "atomic" (40%). The remaining 46% of the sky is 
excluded by our FIRAS mask, or by the lack of Leiden-Dwingeloo survey data at 5 < —30°. 

In Figure [], we plot the correlation slope between our model predictions and the 
FIRAS data for the "molecular" and "atomic" sky. Since the full-sky fits are dominated 
by the "molecular" sky, the model very nearly fits that zone with a correlation slope near 
unity at all frequencies (Figure ||6). However, "the atomic" sky shows deviations relative to 
the model that approach ~ 15% at low frequencies (Figure |c). 

We suspect that these differences between "atomic" and "molecular" zones represent 
true variations in dust grain properties. These variations can be quantified by adjusting our 
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model parameters to achieve a better fit in the "atomic" zone. The fits can be improved by 
retaining the best-fit parameters in the molecular zone, and adjusting a\ lower, f\ higher, 
or <7i/<72 higher in the atomic zone. Lowering a\ does not improve the fit as well, and would 
require a qualitative change in the millimeter opacities of the "silicate" component. We 
consider this possibility the least likely, though it is possible that we are seeing ice mantle 
accretion or some other environment-dependent mechanism. 

In this region of parameter space, f\ and qi/q^ are sufficiently degenerate that an 
almost equally good fit may be obtained in the atomic zone by adjusting either parameter 
a modest amount (Figure f|). Either fi can be increased by ~ 15% or increased by 

~ 25% in the atomic zone (see Figure |6]c). Increasing f\ is equivalent to increasing the 
amount of the "silicate" component relative to the warmer "carbon-based" component. The 
change in q\/q2 could very plausibly be interpreted as a decrease in k\ (an increase in q-i) in 
the molecular zone. If the second component is physically composed of carbon-based grains, 
it might be responsible for the UV absorption bump at 2175 A and sensitive to saturation 
of that feature. Much of the molecular zone has an |SFD98| extinction A(V) ^0.5 mag 
which corresponds to v4(2175 A) ^ 1 mag. When this feature begins to saturate, the carbon 
grains see a change in the ISRF that effectively reduces k\. The small differences we see 
between the "atomic" and "molecular" zones are certainly consistent with spatial variations 
in the radiation field, which we consider to be one of the most reasonable explanations. 
Note that we have not proved that this mechanism is responsible, or even that the two 
dominant components are silicate and carbon-dominated grains. However the data are 
consistent with this interpretation. 

6. APPLICATION TO CMBR ANISOTROPY 

We compare predictions from our best-fit two-component model to the COBE DMR 
data. Significant microwave emission from dust was found by Kogut et al. 1996 by using 
DIRBE 140/iin flux as a dust template. As a similar comparison, we compute the correlation 
slope in each DMR channel with respect to our model. This slope is sensitive only to 
emission correlated with the dust, and does not depend on isotropic backgrounds. For the 
purposes of these fits, DMR pixel i is weighted by the inverse of a { 2 + o cmb where crj is the 
measurement noise in pixel i and ocmb is the RMS power in the CMB anisotropy, taken to 
be 30/iK ( [Bennett et ai. 1996; ). 

The first column in Table [| is the correlation slope of DMR against our best-fit model 
evaluated at 500 GHz. This frequency is chosen because it is low enough to be on the 
Rayleigh- Jeans side of the Planck function, so the dust spectral index between DMR and 
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our 500 GHz predictions has very little temperature dependence. Also, this frequency is 
high enough that FIRAS obtained high quality data for the dust emission. As Figure ^ 
has demonstrated, our model is well tested at 500 GHz and one may be confident that 
it represents real dust emission on the sky. Therefore, 500 GHz is a sensible reference 
frequency to use in such comparisons. The next column of the table is the RMS power in 
/j,K brightness temperature implied by this correlation slope. Note that we are confined to 
the mask used in the computation of the DMR CMBR anisotropy. A less exclusive mask 
would yield a higher RMS power. The remaining columns of Table |3| are similar, but use 
the 140/im DIRBE map instead of our prediction, for direct comparison with the Kogut et 
al. results. However, as we have already seen in Figure 0, the correlation at 500 GHz is not 
tight, and is probably worse at DMR frequencies. These comparisons in Table |^ are meant 
to indicate the spectral shape of dust-correlated microwave emission. 

It is clear in Table [| that the dust-correlated emission at 31 GHz is larger than at 53 
GHZ - when according to our model it should have fallen by a factor of ~ 5. Comparison 
of DMR data with our model evaluated at the same frequency gives an idea of the amount 
of the excess. Again, correlation slopes are computed, for DMR data as a function of 
model predictions, with a slope of unity corresponding to an accurate model. The results 
are tabulated in Table f|. At 90 GHz, there is 20% more dust-correlated emission than 
predicted. At 53 and 31.5 GHz, this emission is too high by a factor of 2.4 and 20, 
respectively. Similar excess emission is seen at 14.5 and 31 GHz from the data of Leitch 
et al. (1997) and at similar frequencies in the Saskatoon data ([Netterfield et al. 1995| ; 



Netterfield et al. 19971 ; |de Oliveira-Costa et al. 19971 ). 



This excess microwave emission is clearly correlated with the dust, but not due to its 
thermal (vibrational) emission. It has recently been suggested that magnetic dipole emission 
from paramagnetic grains ( praine fc Lazarian 1999| ) or electric dipole emission from rapidly 
spinning dust grains ( praine fc Lazarian 1998b| ) could dominate at these frequencies. Others 
have suggested that dust-correlated free-free emission may be responsible ( |de Uliveira-Costa 



et al. 1998|) . However, Draine and Lazarian argue against this on energetic grounds in 



Praine fc Lazarian (1998a]| . Galactic synchrotron emission is not a favored explanation 



because it is unlikely to be highly correlated with the dust. 



6.1. Templates for CMBR Contamination 

It is critical that CMBR experimentalists compare their observations with valid models 
for the Galactic dust emission. A "template approach" is often carelessly used to compare 
observations with expected contaminants, with the correlation amplitude indicating the 
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level of contamination. For example, 100/xm emission maps (e.g., IRAS or DIRBE) or 21 cm 
maps (e.g., Leiden- Dwingeloo: |Hartmann fc Burton 1996| ) are often used as templates 
for microwave dust emission. These templates ignore well-measured variations in dust 
temperature and variations in the dust/gas ratio. We demonstrate this point by differencing 
the broad-band FIRAS 500 GHz map (Figure 0) with best-fit templates. The residuals 
with respect to the H i template (Figure ||) or the 100/im template (Figure |9|) are noticeably 
worse than with our model prediction (Figure [U]). In addition, the Hi or 100/mi template 
offers no insight as to frequency- dependence of the dust emission. Because lOOmu is so 
sensitive to the temperature, it may be the worst of these at high Galactic latitude. 

This need to use the proper template will grow with future data sets. In the case of the 
DMR data, no adverse effects resulted from the use of the 140/xm template, as can be seen 
in Table |3|. The expected RMS power from dust in the DMR channels is not significantly 
altered by using our model. However, the signal-to-noise ratio of DMR is much less than 
1 per pixel. The measurement noise overwhelms the template errors in this case. In the 
case of S/N~ 1 data, e.g. FIRAS 500 GHz, it is apparent from Figure [| that our template 
is more accurate, and much more readily detected, than 100/zm flux or Hi column. These 
considerations will be even more important with satellites such as MAP (13') because our 
map is the only full-sky well calibrated dust model at high resolution (6'). 

We would encourage CMBR researchers to present measurements of dust-correlated 
microwave emission by using our predictions as a baseline, so that (at least where vibrational 
dust emission dominates) the comparison is free of temperature-dependent biases and 
assumptions about the dust/Hi ratio. This will allow easy comparison of samples from 
various parts of the sky - a comparison which is quite difficult with current dust templates. 

7. SUMMARY AND CONCLUSIONS 

We have demonstrated that the SFD98 100/zm emission map, extrapolated with a 
very simple two-component dust model, is an excellent predictor of the Galactic emission 
as seen by FIRAS at all frequencies. Although the older SFD98 v 2 emission predictions 
are tightly correlated with the FIRAS data, the correlation slope deviates significantly 
from unity at frequencies far from 1250 GHz and 3000 GHz (240 and 100/im) where it is 
constrained by the DIRBE data. The v 2 emissivity assumed by SFD98 produces a reduced 
xl ~ 30 when compared with the FIRAS data. Although this fit is unacceptable, no other 
single-component power law emissivity model improves xt significantly. 

We provide a general multi-component model where each component is described by 



-24- 



an emissivity power law, a, power fraction, /, and a ratio of thermal emission to optical 
opacity, q. Each component is required to be in equilibrium with the ISRF. This couples 
the temperatures of each component, both of which are constrained along each line-of-sight 
by the DIRBE 100/240/im ratio. 

In addition to one-component models evaluated for various emissivity power law 
indices, we evaluate 4 two-component models. They correspond to Pollack et ai, Reach et 
ai, our best fit using Pollack et al. emissivities, and our best overall fit for all four model 
parameters. Our best fit parameters are physically reasonable, given empirical evidence 
found in Agladze et al. (1996). See Table [I] for a summary of our results. 

The data argue very strongly that the dust properties of the ISM are uniform 
over virtually the entire high latitude sky on angular scales greater than 7°. We have 
found marginal evidence for variations in the molecular-dominated zones relative to 
atomic-dominated zones. We tentatively suggest that these variations are due to UV optical 
depth effects within the molecular zones. 

This thermal (vibrational) dust emission model fails to explain dust-correlated 
microwave emission observed by DMR. The 90 GHz emission is in approximate agreement 
with our model, but the 53 and 31 GHz DMR emission is high by factors of 2.4 and 20, 
respectively. This excess emission could result from rapidly spinning dust grains ( Praine & 
[Lazarian 1998b ) or from free-free emission. Whatever the emission mechanism, it must be 
strongly correlated with the thermal (vibrational) dust emission. 

Predictions of our best-fit model for thermal dust emission will be made available on 
the World Wide Web.Q 
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A. FURTHER PROCESSING OF THE FIRAS DATA 

The FIRAS pass 4 data products are released in two sets of frequency bins, 
corresponding to the two sides of the instrument. There are 43 bins in the low frequency 
(LOWF) set, running from 68.02 to 639.37 GHz, and 170 bins in the high frequency (HIGH) 
set, running from 612.19 to 2911.29 GHz. The bin spacing is 13.6041 GHz, giving a 3 bin 
overlap between LOWF and HIGH. 



A.l. Recalibration 



An error in the FIRAS external calibrators is described in [Mather et al. 1999| . The 
thermometers were found to be miscalibrated by 5 mK, causing a systematic error in 
LOWF. It is expected that this miscalibration introduces an error of less than one percent 
in the LOWF data, which is negligible for the purposes of this paper. Therefore we have 
chosen to ignore this problem in the data. 

Comparisons of the high-frequency FIRAS data to the DIRBE 240/im data shows an 
inconsistency at the 1% level. The DIRBE gain is uncertain at this level owing to the 
uncertainties in its filter response, and calibration technique. However, we have chosen to 
reduce all the FIRAS measurements by 1%, which is well within the gain uncertainty of the 
HIGH data, and within the measurement noise of the LOWF data. Because the covariance 
of neighboring FIRAS frequency bins is embodied in the FIRAS covariance matrix, our 
results are only weakly dependent upon this 1% recalibration. 



A.2. Bad Bins 

Line emission from CO, [C I], [N II], [C II], [O I], and CH was detected. We have 
excluded those bins and bins corresponding to O2, H 2 0, and [Si I], even though no emission 
was detected by the FIRAS team. A few bins were excluded because of residuals in the 
mirror transport mechanism (MTM) ghost removal. Several bins are excluded from our 
analysis in the frequency range 639.37 to 680.21 GHz for three reasons: (1) inefficiencies in 
the dichroic splitter, (2) a destructive interference pattern caused by reflection off of the 
plastic holder of one of the optical elements, and (3) an aliased MTM sideband. All of 
these effects taken together overwhelm the signal in these bins and justify their exclusion. 
A summary of the frequency bins excluded in our analyses is found in Table |2|. 



Our analyses make use of 123 frequency bins at 100 < v < 2100 GHz. At lower 
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frequencies, the S/N of the dust emission is less than one. At higher frequencies, the 
absolute error in the FIRAS gain exceeds 2% due to uncertainties in the bolometer 
calibration (see |F1RAS Exp. Supp. 15571 §7.3.2). 



B. COMPUTATIONAL METHODS 

Comparison between the Galactic dust emission as observed by FIRAS and our 
predictions is computationally challenging. The predictions are made at the resolution 
of DIRBE to take full advantage of temperature information on scales of ~ 1°. These 
DIRBE-resolution predictions are then smoothed to the FIRAS beam and compared to 
4376 FIRAS pixels (71% of the sky) at 123 frequencies, for a total of ~ 540, 000 data value 
comparisons. 



B.l. Computation of Spectral Shape 

Aside from the 4 global parameters (fx, qi/q?, a%, a%), the model predictions are made 
only from the DIRBE observations at 100 and 240/im. We make use of the DIRBE 100/zm 



map described in |SFD98| with zodiacal light and the cosmic infrared background removed. 
Our dust temperatures are always derived from a ratio map, 7Z, which is a filtered hoo/hio 
flux ratio in the DIRBE passbands. For a multi-component dust model, this ratio map 
measures the following combination of model parameters: 

^ _ J2k Ki 00 (a k , T pk )I 100 (T pk ) 

J2k ^2io( a k: T p k)l2Ao(Tpk) 

For more than one component, the temperatures are related via equation |TJ. For each 
model, we tabulate 7Z as a function of the warmer component, T 2 , as described by equations 



14] and pl| . We fit a 6-th order polynomial to the curve lnT 2 (ln7?.) for the domain 



10<T 2 <31K. For our best-fit two-component model, 

lnT 2 = 2.872 + 0.2407 lnft + 2.962 x 10" 2 In 2 71 + 4.719 x 10" 3 ln 3 ft (B2) 

+9.320 x 10~ 4 In 4 71 + 1.109 x 10~ 4 In 5 71. 

At each DIRBE pixel p, we read the values of the 100/im flux and the ratio 7Z, which 
in turn recovers Tx and T 2 . The flux at any frequency is then given by equation [l] for a 
one-component model or equation |6] for models with more than one component. 

Conceptually, the DIRBE-based predictions are convolved with the FIRAS beam before 
comparison to the FIRAS data. In practice, such a straight-forward approach proved too 
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computationally expensive when minimizing the residuals over several model parameters. 
The dust temperature variations are rarely large within one FIRAS beam, allowing us to 
make approximations for the temperature distribution within a beam. Let % index the 6144 
FIRAS pixels and j index the DIRBE map pixels. Within each FIRAS beam centered on 
FIRAS pixel i, our predictions can be explicitly expressed as 

F l (u)=J2B ll I,Y(u,n,) (B3) 

3 

where B^ describes the beam pattern (the fractional contribution of pixel j to FIRAS pixel 
i), Ij is the 100/im flux in DIRBE pixel j, and Y(v,lZj) describes the model spectral shape, 

The beam pattern is normalized to unity, 

E % = 1. (B5) 

3 



for each FIRAS pixel. Equation |B3| can be rewritten as 

F i (u) = (I i ) J £W ij Y(u,H j ) (B6) 

3 

where we have defined a weighted mean for the 100/im flux, 

(I i )=J2B ij I j , (B7) 

3 

and a weighting function Wij defined as 

W i3 = (B8) 
This weighting function is also normalized to unity within each beam: 

£W<i = l. (B9) 

3 



A direct evaluation of would work, but is very expensive, so we resort to a Taylor 
expansion. 
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B.2. Taylor Expansion 



At each frequency u, Y depends only upon the 100/im/240/im flux ratio, 1Z. We 
expand Y(1Z) about the weighted mean ratio in FIRAS pixel i, IZi, as follows: 



Y{Kj) = Y{%)+Y\%){K j -K i )+-Y"{K i 



Hi - Hi 



TZa - Ki 



+• 



(BIO) 



where the derivatives are with respect to 1Z and we have dropped the v subscript for clarity. 
Computing the weighted sum of Y within one FIRAS pixel yields 



E WvYiKj) = Y(%) + l -Y"(K^(K) + ^'"(iz^iz) + ■ 



where the term linear in (TZj — IZi) vanishes, and af(lZ), 



(Bli; 



(B12) 



is a weighted variance within FIRAS pixel % and 



(B13) 



is a weighted difference cubed. Combining equations and |B11| the flux at any frequency 
is recovered via 



Fi(u) 



(I<) 



(B14) 



Note that this expansion is implemented to describe temperature fluctuations within a 
FIRAS pixel. The values of (Jj), 7Zi, erf (11) and sf(lZ) need only be computed once for all 
the DIRBE values within each FIRAS pixel. Once these values are saved, there is no need 
to return to the higher-resolution DIRBE maps. The predictions for a given dust model 
establishes the relationship Y(1Z), and the flux in FIRAS pixel i is quickly computed via 
equation B14| . 

We have carried the Taylor expansion to third order to establish convergence: the 
third-order terms are significantly smaller than the second-order terms, and are usually 
negligible. All results in this paper are obtained with the third-order Taylor series. Setting 
<J 2 (1Z) equal to zero would ignore these small-scale temperature variations, and introduce 
errors at the level of a few percent in our model predictions. 
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B.3. Definition of x 2 

The comparison between predictions and the FIRAS data is further simplified by 
collapsing the problem spatially. At each frequency, a regression line is computed for the 
FIRAS flux as a function of the predicted flux (as in Figure |2|). The pixel weights from the 
FIRAS data are used for these regressions. The zero-point of the best-fit slope is ignored, as 
it is sensitive to uncertainties in the zodiacal light model or the cosmic infrared background 
(CIB). The slope of the regression is our measure of goodness-of-fit for a model, with a 
slope of unity at all frequencies corresponding to perfect agreement. 

The x 2 f° r each model is computed from the 123 slope values and significance of their 
deviations from unity. The full covariance matrix CV, QF1RAS Fxp. Supp. 19971 , §7.1.2) 



for the FIRAS data is used to couple the errors between frequency bins. We define a 
dimensionless covariance, 

g« = -fa- (B15) 

where i, j index the 123 used frequency bins. The variance in the correlation slope m 
at frequency i, <r 2 (mj), is derived from the linear regression for each frequency, assuming 
uncorrelated Gaussian measurement noise. These a 2 {mi) values are dimensionless, because 
the rrii are dimensionless slopes of order unity. Because of the frequency covariance, 
the variance at frequency % contains contributions from the measurement errors at all 
frequencies j as: 

2/™ \ _ J2 



a im 



E^KOq,- (Bi6) 



This covariance matrix does not include the contribution from the overall bolometer gain 
errors Jj (termed JCJ errors in [FIRAS Exp. Supp. 1997 , §7.3.2). The full covariance matrix 



includes the JCJ terms JiJj, yielding a x 2 of 

X 2 = mi\C[f{mi)o{mj) + JiJ^rrij. (B17) 
This expression for x 2 is used for the fits in this paper. 



C. DATA PRESENTATION 

We provide an electronic data distribution that computes thermal emission from 
Galactic dust for any of the models considered in this paper. The preferred model 
is the two-component model with a.\ = 1.67, Oi2 = 2.70. Intensities are computed at 
any frequency using equation [i] for single-component models and ^ for two-component 
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models. The sky brightness is computed in units of MJy/sr which can be multiplied by 
4024(z//90 GHz)~ 2 yuK( MJy/sr) -1 to convert to a brightness temperature in fiK. Brightness 
temperature may be converted to thermodynamic AT by multiplying by the "planckcorr" 
factor: 

(e x - I) 2 

Planckcorr = ^ — — -L (CI) 

x 2 e x v ' 

where x = hv/kbTcMB and Tqmb = 2.73 In Table |, these factors are evaluated for a 
number of frequencies typical of CMB anisotropy experiments. 

The predictions for thermal emission discussed in this paper are based upon the DIRBE 
100/im map (with zodiacal light and CIB removed) and a ratio map, TZ, which is a filtered 
hoo/hw flux ratio in the DIRBE passbands. These maps are stored as simple FITS images 
in pairs of 1024 x 1024 pixel Lambert ZEA (Zenithal Equal Area) polar projections, similar 
to the data format used for (3FD98 . The NGP projection covers the northern Galactic 



hemisphere, centered at b = +90°, with latitude running clockwise. The SGP projection 
covers the southern Galactic hemisphere, centered at b = —90°, with latitude running 
counterclockwise. (Note that Figs. 0, ||, |, and [It] show the SGP projections rotated by 
180°.) Galactic coordinates (l,b) are converted to pixel positions (x,y) via 

N i (N-l) 

x = — y/l - n sin (b) cos (/) + V J (C2) 

nN i (N-l) 

y = - — y/l-nsm(b) sin (I) + ^-L (C3) 

where = 1024 and n = +1 for the NGP, and n = — 1 for the SGP. Pixel numbers are 
zero-indexed, with the center of the lower left pixel having position (x,y) = (0,0). These 
Lambert projections are minimally distorted at high Galactic latitudes, with the distortion 
approaching 40% at b = 0°. The pixel area of (9f49) 2 oversamples the FWHM of 40'. 

Predictions can be made at higher resolution by extrapolating from the IRAS rather 



than the DIRBE 100/mi map. We use the high-resolution 100/zm map from |SFD98| , 
which contains reprocessed IRAS/ISSA data recalibrated to DIRBE. These maps contain 
4096 x 4096 pixels (N = 4096). The pixel size of (2:372) 2 well samples the FWHM of 
6'1. This map has ~ 20,000 IRAS sources removed, which is appropriate for microwave 
predictions since IR-luminous stars and galaxies are not expected to contribute significantly 
to the microwave sky brightness. 

The caveats to using these maps to predict emission from Galactic dust can be 
summarized as follows: 

1. At frequencies on the Wien tail of the emission (A ^ 100/xm), we underestimate the 
dust emission by not including the contribution from small, transiently heated grains. 
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2. At microwave frequencies [y ^ 100 GHz), we have ignored magnetic dipole emission 
(c.f., praine fc Lazarian 1999| ) and electric dipole emission from rapidly rotating 
grains (c.f., praine fc Lazarian 1998b| ). Either of these mechanisms may be expected 



to dominate the thermal emission. 

3. Our best-fit model is not a complete description of the dust properties. This 
model shows residuals that correlate with such environmental properties as the dust 
temperature and dust/gas ratio. 

4. Unresolved infrared-luminous Galactic sources (primarily stars) are removed from the 
IRAS maps to a flux level of /ioo ~ 0.3 Jy (see [SFD98|) . These stars are not expected 
to contribute significantly to the sky brightness at frequencies v ^ 1000 GHz, but this 
has not been explicitly tested. 

5. Although the angular resolution is 40' for the DIRBE 100/im map and 6fl for the 
reprocessed IRAS 100/im map, our extrapolations to other frequencies relies on an 1Z 
map with an effective resolution of l.°3. 



The data files and corresponding software will be available in the CD-ROM series of the 
AAS, or from our web site.Q Mask files are also available that contain the most important 
processing steps for any given position on the sky. Further details will be available with the 
data files. 
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Fig. I. — Difference spectra from COBE data, after CMBR monopole and dipole removal. 
The bright and faint regions of the sky are differenced for each channel in the DMR 
(diamonds) and FIRAS (solid lines) data sets, excluding the Galactic plane and Magellenic 
Clouds. The sky is divided into cold, warm and hot zones based upon DIRBE /ioo/-^24o color 
ratios. The differences in each zone are renormalized to a 100/xm flux of 1.0 MJy/sr, which 
is a typical flux level for the high-latitude sky. Note the factor of two difference between the 
cold and hot zones at v 700 GHz, relative to the 100/xm normalization. For comparison, 
the dotted line represents 10~ 5 the level of the CMBR spectrum. 
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Fig. 2. — FIRAS - DIRBE comparison. Comparison of FIRAS emission in a synthesized 
500 GHz broad-band versus (a) Hi emission, (b) DIRBE lOOyum (3000 GHz) emission 
with zodiacal contamination removed, (c) prediction from SFD98 using single-component, 
z/ 2 -emissivity model, and (d) prediction from our best-fit two-component model. The 
comparisons are made over 71% of the sky. Straight lines are fit and overplotted using 
the statistical errors in the FIRAS data. The scatter about this line is ~ 3.5 times smaller 
in (c) or (d) as compared to (a) or (&). The slope in (d) is almost unity, as expected for a 
good prediction. 
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Fig. 3. — Dust-correlated emission, scaled by v 2 B v (T) for ease of comparison. The FIRAS 
data (error bars) would be consistent with unity if the v 2 emissivity model were correct. 
Panel (a) overplots broadened temperature models with AT FWHM = 3 K (dashed line) and 
ATfwhm = 6 K (dash-dot line). Panel (b) overplots single-component models with u 15 
(top), v 1 - 7 , and z/ 2 2 emssivity laws. The horizontal dotted line corresponds to z/ 2 . Panel (c) 
overplots two-component models, with the best-fit model shown as a solid line. See Table [1] 
for the specific model parameters. These results are not sensitive to an isotropic background 
in the FIRAS data. The DMR 90 GHz measurement is shown as a diamond. The DMR 30 
and 53 GHz measurements fall well above any model curves. 
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Fig. 4. — Contour of x 2 f° r parameters fi and qi/qi, fixing the emissivity laws (a±, a 2 ) to 
their best-fit values. Our best-fit two-component model is denoted by an X. 
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Fig. 5. — FIRAS versus best-fit model correlation slopes. The sky is divided into three zones 
based upon temperature: (a) cold regions (1Z = /ioo/xm/^24o/im < 0.62), (b) warm regions 
(0.62 < 1Z < 0.69), and (c) hot regions (1Z > 0.69). The systematic residuals between 
zones is not more than ~ 5%. The vertical line is drawn at 240/im, where the models are 
constrained to fit the DIRBE measurements. 
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Fig. 6. — FIRAS versus best-fit model correlation slopes. The sky is divided into two zones 
based upon dust/gas ratio: (a) both atomic- and molecular-dominated zones, (b) zones 
dominated by molecular clouds, and (c) zones dominated by atomic gas. The maximum 
deviation is ~ 15% at 100 GHz for the atomic zone. For the atomic gas, we overplot our 
best-fit model modified to fi = 0.0465 (solid line) or modified to qxjqi = 15.0 (dashed line). 
The vertical line is drawn at 240/xm, where the models are constrained to fit the DIRBE 
measurements. 
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THIS FIGURE AVAILABLE at |ittp://astro.berkeley.edu/dust/index.html 

Fig. 7. — FIRAS broad-band 500 GHz map, as defined in the text (equation These 
Lambert ZEA polar projections are centered on the NGP (left), and SGP (right), with 
Galactic latitude labeled in degrees. Lines of constant latitude and longitude are spaced 
every 30 degrees. 



THIS FIGURE AVAILABLE at [http://astro.berkeley.edu/dust/index.html 



Fig. 8. — Difference map between the broad-band FIRAS 500 GHz map and the best-fit 
Hi template convolved to the same beam shape. Declinations south of 5 = —30° were not 
observed in the Leiden-Dwingeloo (Hartmann & Burton 1996) survey, accounting for the 
blank regions of missing data. The thin black lines outline FIRAS pixels masked from our 
analysis (see § |2.4.2|) . 



THIS FIGURE AVAILABLE at [http://astro.berkeley.edu/dust/index.html 



Fig. 9.— Difference map between the broad-band FIRAS 500 GHz map and a best-fit 
lOOyum template convolved to the same beam shape. 



THIS FIGURE AVAILABLE at |http: / /astro. berkeley.edu/dust/index.html 

Fig. 10. — Difference map between the broad-band FIRAS 500 GHz map and the emission 
predicted from our best-fit two-component dust model convolved to the same beam shape. 
The residuals are noticeably smaller than residuals obtained by using either an H i template 
(Figure |8|) or 100/rni template (Figure 0). 
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Table 1. Fit results for dust emission models 



# 


Model 




«2 


h 


<7i/<72 




(T 2 ) 


Pi /Pi 


x 2 




1 


One-component: v emis 


1.5 




1.0 


1.0 


20.0 






24943 


204 


2 


One-component: v emis 


1.7 




1.0 


1.0 


19.2 






8935 


73 


3 


One-component: ^ 2 emis 


2.0 




1.0 


1.0 


18.1 






3801 


31 


4 


One-component: v 22 emis 


2.2 




1.0 


1.0 


17.4 






9587 


79 


5 


Pollack et al. 2-component 


1.5 


2.6 


.25 


0.61 


17.0 


17.0 


.33 


1866 


15.3 


6 


Two-component: both v 2 


2.0 


2.0 


.00261 


2480 


4.9 


18.1 


.0026 


1241 


10.3 


7 


Two-component: fit /, q 


1.5 


2.6 


.0309 


11.2 


9.6 


16.4 


.0319 


244 


2.03 


8 


Two-component: fit f,q,a\,a% 


1.67 


2.70 


.0363 


13.0 


9.4 


16.2 


.0377 


219 


1.85 



Note. — The dust models are described by ai, Q2, fi, and q\/qi- The mean temperatures for each dust 
component, (Ti) and (T2), are evaluated for the mean hoo/hio color ratio in the high-latitude sky. The ratio 
of power emitted by each component is Pi /P 2 . 
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Table 2. FIRAS bad frequency channels 



Reason 


Channel 


V 


V 


A 






( cm- 1 ) 


( GHz) 


(H 


CO (J=l-0) 


3 


3.631 


108.83 


2753.8 


CO (J=2-l) 


12 


7.717 


231.27 


1295.9 


MTM 


18 


10.440 


312.90 


957.8 


MTM 


19 


10.894 


326.50 


917.9 


CO (J=3-2) 


20 


11.348 


340.10 


881.2 


o 2 


26 


14.072 


421.73 


710.6 


CO (J=4-3) 


29 


15.433 


462.54 


647.9 


[CI] 


31 


16.341 


489.75 


611.9 


H 2 


36 


18.611 


557.77 


537.3 


CO (J=5-4) 


37 


19.065 


571.37 


524.5 


Dichroic 


40 


20.427 


612.19 


489.6 


llph T*f~ll C 


41 


90 881 

ZjU .go ± 


695 79 


478 9 


Dichroic 


42 


21.335 


639.40 


468.7 


Dichroic 


43-48 


21.6 


646 


464 


CO (J=6-5) 


49 


23.150 


693.81 


432.0 


CO (J=7-6) 


57 


26.782 


802.65 


373.4 


[CI] 


58 


27.236 


816.25 


367.2 


H 2 


80 


37.222 


1115.54 


268.7 


[Nil] 


105 


48.570 


1455.65 


205.9 


H 2 (J=2-l) 


124 


57.195 


1714.12 


174.8 


[CII] 


137 


63.096 


1890.98 


158.5 


[CII] 


138 


63.550 


1904.58 


157.4 


[CII] 


139 


64.004 


1918.19 


156.2 


[CII] 


140 


64.457 


1931.79 


155.1 


[01] 


149 


68.543 


2054.23 


145.9 


[Si I] 


168 


77.167 


2312.71 


129.6 


[Nil] 


179 


82.161 


2462.35 


121.7 


CH (J=2-l) 


188 


86.246 


2584.79 


115.9 
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Table 3. Comparison of dust templates with DMR 



V 


Model(500 GHz) 


RMS 


DMR/140/zm 


RMS 


GHz 


slopexl0~ 3 




slopexl0~ 4 


fiK 


31.5 


1.81 ±0.27 


20.5 ±3.1 


1.81 ±0.28 


20.6 ±3.2 


53 


1.35 ±0.29 


5.4 ±1.1 


0.99 ±0.29 


4.0 ± 1.2 


90 


4.27 ± 1.08 


5.9 ± 1.5 


3.38 ± 1.10 


4.7 ± 1.5 



Note. - Correlation slopes of the DMR channels with best- 
fit 500 GHz prediction and 140^m DIRBE data. Slope values are 
dimcnsionlcss flux ratios. RMS values are the RMS power, in /xK, 
expected in the DMR maps due to the dust emission traced by 
template. These values compare to an RMS power of 29 ± l^K due to 
CMBR anisotropy. 
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Table 4. Excess dust-correlated DMR emission 



V 


DMR/Model(z/) 


% variance 


RMS power 


GHz 








31.5 


20.09 ±3.09 


0.87 


20.2 ±3.1 


53 


2.45 ±0.50 


0.57 


5.6 ±1.1 


90 


1.18 ±0.29 


0.42 


6.2 ± 1.5 



Note. - Excess dust-correlated microwave emission 
measured by DMR. 1) frequency of DMR channel, in 
GHz. 2) correlation slope of DMR emission vs. model 
predictions. 3) percent of variance in DMR data 
accounted for by this dust. Note that the vast majority of 
the variance is receiver noise. 4) RMS power due to dust 
emission, in fiK brightness temperature. 
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Table 5. DIRBE K-correction fit coefficients 



coefficient 


a = 1.50 


a = 1.67 


a = 1.70 


a = 2.00 


a = 2.20 


a = 2.60 


a = 2.70 


a 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


1.00000 


ai 


2.08243 


2.15146 


2.14106 


2.18053 


2.55941 


3.16383 


3.31600 


a-2 


-4.72422 


-4.84539 


-4.83639 


-4.89849 


-5.41290 


-6.23131 


-6.43306 


a-3 


2.29118 


2.35210 


2.35919 


2.38060 


2.57867 


2.86900 


2.93939 


bo 


-0.88339 


-0.87985 


-0.93625 


-0.80409 


-0.80318 


-0.50356 


-0.41568 


h 


4.10104 


4.10909 


4.19278 


3.95436 


4.20361 


4.07226 


4.02002 


b 2 


-4.43324 


-4.43404 


-4.46069 


-4.27972 


-4.55598 


-4.70080 


-4.72432 


b 3 


1.76240 


1.76591 


1.77103 


1.70919 


1.80207 


1.87416 


1.88865 



Note. — Fit coefficients for DIRBE 100/im band color correction factors Kioo(a, T), fitted using 
equation @. 
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Table 6. Unit Conversion Factors for Selected CMBR Experiments 



Experiment 


V 


A 


factor 


Planckcorr 




GHz 


mm 


/zK/(MJy/sr) 




COBE/DMR 


31.5 


9.52 


32849 


1.02582 




53.0 


5.66 


11603 


1.07448 




90.0 


3.33 


4024 


1.22684 


MAP 


22.0 


13.64 


67344 


1.01253 




30.0 


10.00 


36216 


1.02340 




40.0 


7.50 


20371 


1.04190 




60.0 


5.00 


9054 


1.09623 




90.0 


3.33 


4024 


1.22684 


Planck 


30.0 


10.00 


36216 


1.02340 




44.0 


6.82 


16836 


1.05087 




70.0 


4.29 


6652 


1.13275 




100.0 


3.00 


3259 


1.28562 




143.0 


2.10 


1594 


1.65110 




217.0 


1.38 


692.2 


2.98186 




353.0 


0.85 


261.6 


12.8186 




545.0 


0.55 


109.7 


157.85 




857.0 


0.35 


44.38 


15392 



Note. Column 1 contains the frequency, in GHz, for 

which the unit conversion factors are computed. Column 2 
is the corresponding wavelength in mm. A value in units of 
MJy/sr should be multiplied by the factor in Column 3 to 
convert to fiK brightness temperature. Brightness temperature 
is multiplied by Planckcorr (Column 4; equation |ClD to convert 
to thermodynamic temperature, assuming Tcmb = 2.73K. 



